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Abstract. We investigate the local density of states of the one-dimensional half- 
filled spinless fermion model with nearest-neighbor hopping t > and interaction V 
in its Luttinger liquid phase — 2i < V < 2t. The bulk density of states and the local 
density of states in open chains arc calculated over the full band width ~ At with an 
energy resolution < O.OSi using the dynamical density-matrix renormalization group 
(DDMRG) method. We also perform DDMRG simulations with a resolution of O.Oli 
around the Fermi energy to reveal the power-law behaviour -D(e) |e — epj" predicted 
by the Luttinger liquid theory for bulk and boundary density of states. The exponents 
a are determined using a finite-size scaling analysis of DDMRG data for lattices with 
up to 3200 sites. The results agree with the exact exponents given by the Luttinger 
liquid theory combined with the Bethc Ansatz solution. The crossover from boundary 
to bulk density of states is analyzed. We have found that boundary effects can be seen 
in the local density of states at all energies even far away from the chain edges. 
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1. Introduction 

One-dimensional conductors have fascinated physicists for more than 50 years [U [2] 
because they feature unusual properties which set them apart from ordinary metals. 
Our understanding of ordinary metals is based on the Fermi liquid paradigm [3]. In one 
dimension, however, this theory fails. Instead, the low-energy physics of one-dimensional 
conductors is described by the Luttinger liquid paradigm [21 HI [5]. The predictions of 
the Luttinger liquid theory differ fundamentally from those of the Fermi liquid theory. 
For instance, Fermi liquids have a finite density of states D{e) at the Fermi energy ep. 
In contrast the bulk density of states of Luttinger liquids vanishes as a power law at the 
Fermi energy 

i:>(e) ~ le-epl" (1) 

where the exponent a > depends on the system. This feature is regarded as one 
hallmark of a Luttinger liquid. 

The density of states D{e) can be measured experimentally in photoemission 
spectroscopy and in scanning tunnelling spectroscopy (STS). The STS method yields a 
spatially-resolved local density of states (LDOS) with a resolution of a few A. The 
experimental observation of vanishing densities of states has been reported in the 
photoemission or STS spectrum of various quasi-one-dimensional conductors such as 
Bechgaard salts [6], the organic charge transfer salt TTF-TCNQ large samples 
of single- walled carbon nanotubes [Hj, and the purple bronze Lio.gMoeOiy [HI ITO] . 
Consequently, these materials are believed to be realizations of Luttinger liquids 
although this interpretation remains often controversial. Only very recently, the power- 
law behaviour ([1]) has been observed unambiguously in the STS and photoemission 
spectra of gold wires deposited on semiconducting Ge(OOl) surfaces [TTj . 

From a field-theoretical point of view, the low-energy properties of Luttinger 
liquids are very well understood thanks to powerful methods such as bosonization 
and renormalization group [21 Hj [5]. However, field theory only describes the low- 
energy scaling |e — ep] — )■ of various physical quantities such as the density of states. 
Experimentally, this asymptotic behaviour could only be observed in a more or less 
broad window of excitation energies. Indeed, as real materials are three-dimensional, 
there is always a dimensional crossover [H [2] at low excitation energy below which 
one-dimensional physics can no longer be observed. Moreover, as field-theoretical 
investigations are based on the linear dispersion of excitations close to the Fermi energy, 
there is always a high-energy limit above which the density of states should deviate 
from the power-law ([1]) because of the band curvature in real materials [12]. Therefore, 
the theory must be extended to finite-energy scales beyond the asymptotic behaviour 
covered by the Luttinger liquid theory to facilitate the interpretation of experiments in 
one- dimensional conductors. 

We know that the Luttinger liquid paradigm describes the low-energy properties 
of various one- dimensional quantum lattice models for interacting electrons in their 
metallic phases, such as the Hubbard model [13] away from half filling and the 
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spinless fermion model [21 Lattice models enable us to study the influence of 
finite energy scales such as the curvature and finite width of excitation bands in one- 
dimensional conductors. Some of these models are exactly solvable by the Bethe Ansatz 
method. However, it is very difficult to obtain dynamical correlation functions related 
to spectroscopic experiments from a Bethe Ansatz solution. (The spectral functions 
of the spinless fermion models have been investigated only very recently using the 
Bethe Ansatz [H] but the issue of the density of states has not been discussed in that 
work.) Therefore, the spectral properties of these models at finite excitation energy 
have mostly been determined using numerical methods for the quantum-many body 
problem [15] such as exact diagonalizations, quantum Monte Carlo simulations and the 
density-matrix renormalization group (DMRG) method. Several power-law divergences 
(~ |e|" with a < 0) predicted by field theory have been confirmed with these numerical 
methods. Yet the power-law singularity (II]) with an exponent a > has not been 
observed univocally at finite-energy scale in quantum lattice models so far. 

The DMRG method is one of the most powerful numerical method for computing 
the properties of one- dimensional quantum lattice models [TS] [16]. The density 
of states of some lattice models have been investigated using the original DMRG 
method [TTJ [HI [19] . In these studies DMRG was used to compute the spectral weight of 
the lowest few eigenstates in order to verify the prediction of the bosonization approach 
in the asymptotic low-energy limit. The dynamical DMRG (DDMRG) is an extension 
of DMRG which makes possible the calculation of dynamical correlation functions over 
their full band width [201 [21] . It yields spectra which are broadened by a Lorentzian of 
width 7] which sets the actual energy resolution. Over the last decade DDMRG has been 
used successfully in many studies of spectral properties in quantum lattice models with 
energy resolution down to a few hundredths of the bare band width |22]. Nevertheless, 
this accuracy has not allowed for a direct observation of the power-law behaviour at 
finite excitation energy until now. For instance, DDMRG has been used to calculate the 
complete single-particle spectral functions of the Hubbard model away from half filling 
with a resolution of ?7 = O.lt [22]. The dispersion of holon and spinon branches in the 
DDMRG spectral functions agree perfectly with the exact Bethe Ansatz dispersions. Yet 
at the Fermi energy, where the power law ([l]) should be seen, the momentum-integrated 
spectral weight barely shows a shallow dip. 

In this paper we investigate the density of states of the half-filled spinless fermion 
model in its Luttinger liquid phase. In this model we can study significantly larger 
systems and thus reach a much better resolution than in the Hubbard model. Moreover, 
in the half-filled spinless model any exponent a > can be achieved while in the Hubbard 
model only 0<a;<l/8is possible. 

Both the bulk density of states and the local density of states close to a chain 
edge are calculated numerically using the DDMRG method. They are determined over 
the full band width ~ 4t with a resolution of 77 = 0.04t or O.OSt. Using lattices with 
up to 1600 sites we have also performed DDMRG simulations with a resolution of 
f] = O.Olt around the Fermi energy to reveal the power-law behaviour ([T|) predicted by 
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the Luttinger liquid theory for bulk and boundary density of states. The exponents a 
are determined using a finite-size scaling analysis of DDMRG data for lattices with up 
to 3200 sites and compared to the predictions of the Luttinger liquid theory combined 
with the exact Bethe Ansatz solution. Finally, we discuss the crossover from boundary 
to bulk density of states as one moves away from the chain edges. 

Our paper is organized as follows: In the next section we introduce the model and 
method used in this work. In the third section we discuss our results for the bulk density 
of states while the LDOS close to a chain edge is analyzed in the fourth section. Finally, 
our findings are summarized in the last section. 



2. Models and method 



The one-dimensional spinless fermion model is one of the simplest realizations of a one- 
component Luttinger liquid in a lattice model. It can be interpreted as a system of 
spin-polarized electrons. The model is defined by the Hamiltonian 

N-l JV-l . 1\ / 1\ 

^ = E + sVis) + ^ E (^^- - 2) y"^^' - 2) 

where c| and c- are the creation and annihilation operators for a spinless fermion on 
site j, and the density operator is nj = CjCy The model parameters are the hopping 
amplitude t > between nearest-neighbor sites and the Coulomb interaction V between 
particles on nearest-neighbor sites. The half- filled system corresponds to N/2 fermions 
on the A^-site lattice. 

The one-dimensional spinless fermion model is exactly solvable by the Bethe Ansatz 
method [21 [5]. For —2t < V < 2t its excitation spectrum is gapless and its low-energy 
properties are described by the Luttinger liquid theory. At half filling the dispersion of 
elementary excitations [23] is given by 

e{k) = 2t*\sm{ka)\ (3) 

with the renormalized "hopping term" 



r - f (4) 

2 arccos [^^) 

The universal properties of a one-component Luttinger liquid are determined by two 
parameters: The velocity of elementary excitations (renormalized Fermi velocity) v and 
a dimensionless Luttinger parameter K. From the Bethe Ansatz solution we know the 
relation between these Luttinger liquid parameters and the lattice model parameters at 
half filling, 

2at* 

v = ^ (5) 

with t* given by and 

TT 1 . , 
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where a is the lattice constant. In the non-interacting chain = this yields K = 1 and 
the usual Fermi velocity v = vp = 2ta/h of the one-dimensional tight-binding model. 
The local density of states (LDOS) is defined by 

2 

5{e - En + Eq) for e > ep 
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5{e + En - Eq) for e < ep 



(7) 



where \n) denotes the eigenstates of the Hamiltonian H and E^ their eigenenergies in 
the Fock space. The ground state for the chosen number of particles correspond to 
n = 0. The total spectral weight is 

deD{j,e) = l. (8) 

As the Hamiltonian ([2]) is invariant under the particle-hole transformation o Cj (— !)•' , 
its density of states is symmetric at half filling: e) = D[j, — e) with ep = 0. 

We use the DDMRG method [201 EI] to compute the density of states DDMRG 
simulations yield the convolution of a Lorentzian of width 77 > with the local density 
of states on a A^-site lattice 

As we are interested in the density of states in the thermodynamic limit, we should, 
in principle, extrapolate our numerical data to an infinite lattice size A^ — )■ 00 and 
then to a vanishing broadening 77 — )■ 0. However, the spectrum of a A^-site chain 
is indistinguishable from the same spectrum in the thermodynamic limit if they are 
compared on an energy scale (resolution) Ae > W/N . The parameter W has to be 
determined empirically but it is typically of the order of the effective band width of 
the spectrum considered (here W ~ 4t*). Thus if we choose a broadening 77 > W/N, 
DDMRG spectra can be regarded as the spectra of infinite systems with the same 
broadening 77, which sets the actual resolution. In this work we have systematically 
used 

The DMRG algorithm truncates the Hilbert space in an optimal way to represent 
some chosen target states such as the ground state. Thus DMRG data are affected by 
truncation errors which depend on the number of density-matrix eigenstates kept in the 
calculation [151 CH] • In the DDMRG method the truncated Hilbert space is optimized 
to represent the ground state and eigenstates with a given excitation energy e or lying 
in a given excitation energy range [ei, 62] with £2 — ^ The calculations are carried 
out independently for each excitation energy. We have kept up to 400 density-matrix 
eigenstates in our simulations. Numerical errors are always very small in absolute values 
or when they are compared to the total spectral weight ([8]). Thus the DDMRG results 
presented here can be regarded as numerically exact. This is illustrated in the following 
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Figure 1. Bulk density of states for V = —l.9t { 

( ) calculated in the middle of a 400-site chain (77 



), V^ = 0(- 
= 0.04t). 



■), and V = 2t 



sections by comparison with exact spectra for non-interacting fermions. However, as the 
spectral weight approaches zero at the Fermi energy, relative errors become unavoidably 
larger and can be seen as a scattering of DDMRG data close to the Fermi energy. In 
our numerical calculations we use the energy scale t = 1, the time scale h/t = 1 and the 
length scale a = 1. 

3. Bulk density of states 

In principle, the bulk density of states D{e) should be calculated as the average of the 
LDOS ([7]) over the full lattice in the thermodynamic limit. Here we identify the bulk 
density of states with the LDOS in the middle of a A^-site chain, D{j = N/2, e). While 
this seems intuitively correct for insulators, in a Luttinger liquid boundary effects could 
a priori be felt as far as the chain center. First, we have verified the validity of this 
approach for several exactly solvable non-interacting models. For instance, the bulk 
density of states in a non-interacting chain (V^ = 0) can be calculated analytically. 



for |e| < 2t and D{e) = for |e| > 2t. Figure [T] shows the LDOS calculated with DDMRG 
in the middle of a 400-site non-interacting chain with 77 = 0.04t. On the scale of this 
figure the DDMRG spectrum is indistinguishable from the exact result ffTTl) convolved 
with a Lorentzian of the same width 77. This demonstrates the accuracy of the DDMRG 
method and also confirms that the scaling of the broadening ( fTOi) is appropriate for the 
density of states of the spinless fermion model. 

Additionally, we have checked in interacting systems with DDMRG that D{j,e) 
depends negligibly on the site position j close to the chain center and that D{j = N/2, e) 
is indistinguishable from the momentum-integrated spectral functions [21]. Finally, we 
will examine the crossover from boundary LDOS to bulk LDOS in more detail in the 
next section. All our results confirm that the DDMRG spectra calculated in the middle 
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Figure 2. Position ep/2t of the peaks in the bulk density of states calculated with 

DDMRG in 200-site chains (O ) and renormalized "hopping term" t* /t from (j4| ( ) 

as a function of the interaction parameter V. 



of a spinless fermion chain with a broadening (ITO!) can be regarded as its bulk density 
of states in the Luttinger liquid phase. 

Figure [T] shows the bulk density of states over the full band width for three different 
values of the interaction V. The band edge singularities of the non-interacting density 
of states ( fTTj) are seen as two peaks at e = ±2t because of the broadening rj. Similar 
peaks are visible at positions e = ±ep for V ^ 0. The energy ep increases monotonically 
with V starting from for F = —2t. Actually, figure [2] shows that €p/2t varies exactly 
as the effective "hopping term" t* (jlj) as a function of V. Thus ep = 2t* and the peak 
positions correspond to the edges of the elementary excitation bands [i.e., the maximum 
of the dispersion ([3])] for all values of V. 

In figure [1] we see some spectral weight beyond ep which is clearly due to the 
broadening t] for ^ = but cannot be explained by this broadening for V = —1.9t 
and V = 2t. Therefore, excitations which are made up from more than one elementary 
excitations contribute to the density of states ([7]) at energy |e| > ep. 

According to the Luttinger liquid paradigm the bulk DOS should vanish as the 
power law at the Fermi energy. The Luttinger liquid theory also predicts that the 
exponent is 

for a one-component Luttinger liquid. Thus for the spinless fermion model this exponent 
can be calculated exactly as a function of the interaction parameter V using the Bethe 
Ansatz relation for the Luttinger parameter ([6]). In the non-interacting chain [V = 0) 
this yields a = 0, which implies that the bulk DOS does not vanish but approaches a 
constant value at the Fermi energy, while one gets a > for interacting chains (V ^ 0). 

The power-law behaviour is not visible on the scale of figure [H The vanishing of 
the density of states is just suggested by narrow minima which are seen at e = for 
interacting fermion systems. Figure [3] shows the bulk density of states around the Fermi 
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Figure 3. Bulk density of states around tlie Fermi energy for V ~ 2t calculated with 
DDMRG in a 1600-site chain (77 = O.Oli). The soUd line is the power law ^ with the 
exact exponent a = j. 

energy at a higher resolution (77 = O.Olt) for V = 2t. The power law ([1]) is also shown 
with the exact exponent a = ^ for \^ = 2t. Clearly, the DDMRG data agree with the 
power law over the energy range shown in this figure. 

Actually, we have often found an approximate power-law behaviour over a rather 
broad energy range of the order of t* when the exponent a is not too large. This is 
illustrated in figure H] which shows the bulk density of states for two values of V on a 
double logarithmic scale. In both cases the DDMRG data clearly deviate from the power 
law at high energy e ^ t*. This is not surprising as the curvature of the dispersion ([3]) 
becomes relevant above this energy scale. The DDMRG data can also deviate from 
the power law at low energy € < t] = O.Olt, where the spectrum is dominated by the 
Lorentzian broadening iQ, although this is not apparent in both examples in figure HJ 
For V = 2t the DDMRG spectrum agree particularly well with the power law. For most 




e/t 

Figure 4. Bulk density of states for V = -\f2t (O ) and V = 2t (□) calculated with 
DDMRG in a 1600-site chain (r/ = O.Oli). The solid lines show the power law ([1]) with 
the exact exponent a = j for both cases. 
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Figure 5. Bulk density of states around the Fermi energy for V = — 1.9t calculated 
with DDMRG in a 1600-site chain (77 = O.Olt). The sohd line is the power law ^ with 
the exact exponent a « 1.574 but shifted vertically by 0.0145. 

values of V, however, there is only a rather rough agreement as shown by the second 
example, V = — -\/2t. 

For large exponents a > 1 discrepancies become apparently significant. This is 
illustrated in figure [5] which shows the density of states for V = —1.9t corresponding to 
a ~ 1.574. Clearly, a simple power law ([T]) can not describe the DDMRG data because 
they converge to a finite value for e — )■ 0. This is due to the convolution iQ which 
transfers a spectral weight -DDDMRG(eF = 0) ~ 77 ~ where the original spectrum 
has none. While this finite-size effect occurs for all values of V, it becomes relevant for 
large exponents a only. Nevertheless, a simple shift is enough to offset the finite density 
of states at the Fermi energy and to recover an excellent agreement with the low-energy 
power-law behaviour as shown in figure O 

We note that a similar effect is seen in the STS spectra of atomic chains [H] 
and the purple bronze [9]. Clearly, the experimental densities of states shown in these 
publications are well fitted by power laws around the Fermi energy but they remain 
finite at the Fermi energy. There, interchain couplings, finite temperature or disorder 
may be responsible for this deviation from a pure power-law behaviour. 

Figure H] suggests that the bulk density of states decreases approximately as a 
power law over a rather broad energy range. Yet we have rarely succeeded in extracting 
accurate values for the exponent a from fits to such spectra. The exponent a can be 
determined more accurately and with less computational effort using a finite-size scaling 
analysis [20]. For this purpose we calculate the density of states -D(e) for several system 
sizes N at an energy e which scales as We know that D{e ~ 1/-^) scales as 

A^~° in the thermodynamic limit and thus we can determine a using a power-law fit. 
In this work we have used e = 77 = 16t/N and up to = 3200 sites for this scaling 
analysis. The exponents that we have obtained in the range 0.08 < a < 1.6 agree 
within 10% with the predictions of the Luttinger liquid theory ( 1T2|) combined with the 
Bethe Ansatz solution ([6]). Larger exponents than a 1.6 can be achieved in the half- 
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filled spinless fermion chain for —1.9t > V > —2t. We have not explore this regime 
because the effective energy scale (jlj) becomes rapidly very small and thus numerical 
simulations become increasingly difficult. A reliable determination of smaller exponents 
than a ~ 0.08 would require DDMRG computations for a larger range of system sizes 
than in this work (i.e., N ^ 10'^). In principle, a similar approach could be used to 
determine whether the peaks found at e = ±2t* are band edge singularities like in the 
non-interacting system (|TT|) or smoother structures but we have not carried out this 
analysis yet. 

4. Local density of states 

We now examine the LDOS D{j,e) close to a chain end (hard wall boundary). This 
configuration can also be realized experimentally, for instance in gold chains on Ge 
surfaces [H]. Previous theoretical works [21 [171 UHl EH ES] have shown that in a 
Luttinger liquid the LDOS at chain edges differs significantly from the bulk density 
of states. In particular, the low-energy behaviour is very different for spinless fermions 
with attractive {V < 0) or repulsive {V > 0) interactions. Additionally, the LDOS 
of the spinless fermion model next to a site impurity has been investigated using the 
functional renormalization group (fRG) method [26]. (In a Luttinger liquid the LDOS 
next to an impurity, however weak, is similar to the LDOS close to a chain end in the 
asymptotic low-energy limit [2ll5].) 

We first explore the boundary LDOS, i.e D{j = l,e) on the first lattice site. In a 
semi-infinite non-interacting chain {V = 0) this LDOS is given by 

D{j = l,e) = ^V^t^. (13) 

The DDMRG spectrum for = is shown in figure E] with a broadening t] = 0.04t. On 
the scale of this figure there is no visible difference between these DDMRG data and 
the exact result ( |T3l) with the same broadening. This confirms that the scaling of the 




Figure 6. Boundary LDOS D{j = 1, e) calculated with DDMRG on the first site of a 
400-site chain (77 = 0.04t) for ^ = ( =lt { ), and V ^ 2t ( ). 
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Figure 7. Boundary LDOS D{j — 1, e) calculated with DDMRG on the first site of a 
200-site chain (77 = O.OSt) for F = ( ), V ^ ~lt { ), and V = -IM ( ). 



broadening (fTOj) is appropriate for the LDOS close to chain edges too and demonstrates 
again the accuracy of the DDMRG method. 

Figure [6] also shows the boundary LDOS calculated with DDMRG for two repulsive 
interactions V. The semi-elliptic spectrum ( IT3l) of the non-interacting chain seems to 
split into two smaller semi-elliptic bands which move apart as V increases. A narrow 
minimum (but no gap up to = 2t) appears at the Fermi energy e = 0. The onsets of the 
upper and lower bands correspond to the peak positions e = ±2t* in the bulk density of 
states discussed in the previous section and thus to the edges of the elementary excitation 
bands ([3]). As in the bulk density of states we observe some spectral weight beyond 2t* 
which can be explained by the convolution ([9]) for V = but not for interacting fermions 

In contrast for attractive interactions < the boundary LDOS develops an 
increasingly sharp peak at the Fermi energy as V decreases. Figure [7] shows the DDMRG 
spectra for two values of < as well a.sV = for comparison. For = — t we can still 
see that the spectrum has apparent onsets at energies e ~ ±2t* corresponding to the 
band edges of elementary excitations ([3]). For larger V we can no longer distinguish the 
spectrum onsets but most of the spectral weight is still concentrated within the band of 
elementary excitations —2t* < e < 2t*. 

The Luttinger liquid theory predicts that the boundary density of states also follows 
a power law but the boundary exponent is 

«B = - 1 (14) 

for a one-component Luttinger liquid [21 [17]. Previous DMRG investigations have 
confirmed this expression for the boundary exponent in terms of the Luttinger 
parameter (|6]) obtained from the Bethe Ansatz solution [171 [H] • 

For repulsive interactions > the Luttinger parameter is -R' < 1 according to 
and thus ob > a > 0. This increase of the exponent at a chain boundary has been 
observed experimentally in gold chains on semiconducting Ge surfaces [11] . The power- 
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law behaviour is not visible on the scale of figure Eland the narrow minima at the Fermi 
energy just hints at the presence of a pseudogap. As for the bulk density of states we 
have examined the low-energy LDOS with resolution of 77 = O.Olt and obtained similar 
results. The boundary LDOS follows an approximate power law over an energy range 
~ t* for small exponents a. For a ^ 1 (<(=^ V ~ 2t) we have to offset some finite 
spectral weight at ep = to recover the correct low-energy behaviour. Finally, using 
a finite-size scaling analysis of our DDMRG data we obtain exponents a which agree 
within 10% with the exact results of the Luttinger liquid theory ([14]) combined with the 
Bethe Ansatz solution (Q. 

In contrast for attractive interactions V < the Luttinger parameter is K > 1 and 
thus c^B < 0. Therefore, the Luttinger liquid theory predicts a power-law divergence in 
the boundary LDOS at the Fermi energy. As the effective band width ~ 4t* and the 
exponent ae decrease rapidly with decreasing V, this singularity should become stronger 
and sharper. When V approaches —2t, the peak width must vanish according to (jl]) 
while the exponent ae converges to -1 according to (El) and (HM . Thus the power- law 
singularity turns into a Dirac (5-peak at the Fermi energy in that limit. The predictions 
of the Luttinger liquid theory agree with the boundary LDOS calculated with DDMRG 
which are shown in figure [71 In the DDMRG spectra the Fermi-energy singularity has 
been smoothed into a peak by the broadening (Q. Again we can perform a finite-size 
scaling analysis to confirm that the peak is a singularity and to determine the precise 
values of the exponent a^- This analysis is much easier for power- law singularities with 
negative exponents (i.e., divergences) than with positive ones (pseudogaps) . Thus we 
find exponents which agree very well (within 2%) with the exact values given by the 
Luttinger liquid exponent (1141) combined with the Bethe Ansatz solution ([6]). 

We now turn to the crossover from the boundary LDOS to the bulk density of 
states. First, it is helpful to discuss the LDOS close to the chain edge in a semi-infinite 
non-interacting chain (V^ = 0). Numbering the sites j = 1,2, . . . from the chain edge we 
obtain for small j <^ N and |e| < 2t 



where Tj{x) = cos(arccos(x)j) are the Chebyshev polynomials. As for the bulk density 
of states D{j, e) = for |e| > 2t. The boundary LDOS ( [T3l) is recovered for j = 1. The 
LDOS D[j,e) oscillates between and twice its bulk value ( [TTl) . There are exactly j 
peaks in D{j, e) considered as a function of e. In particular, there is a peak at the Fermi 
energy for odd j but D{j, e) vanishes as for even j. At low energy these oscillations 
correspond to two modes (—1)"' and cos(ej'/t) = cos[2eja/(^f^)]. Field theoretical 
studies have shown that these same two oscillation modes also dominate the low-energy 
LDOS of Luttinger liquids with the renormalized velocity v substituted for vp [Mj . 

As f[T5l) has been derived for j -C A^, it is not surprising that it does not converge 
toward the bulk density of states ( [TTl) for large j. Nevertheless, if we broaden each 
LDOS D{j, e) with a Lorentzian (Q (or a Gaussian, . . . ) of width ~ 4t/j, we find that 
the result converges toward the bulk density of states for j — )■ 00. As in the middle 




(15) 
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Figure 8. LDOS D{j, e) calculated on a 200-site chain [r] = O.OSt) for V = 2t. (a) 
From one edge to the middle of the chain: j — 1, 10, 20, . . . , 100 from top to bottom, 
(b) Close to the chain edge: j — 1,2, ... ,10 from top to bottom. 

of a A^-site chain this broadening ~ At/j is smaller than the one used in DDMRG 
calculations ffTOj) . this property explains why the LDOS calculated in the middle of a 
chain with DDMRG agrees so well with the (broadened) bulk density of states (fTTl) . 

The crossover from boundary to bulk density of states is shown in figure [8] for a 
repulsive interaction V = 2t and in figure [9] for an attractive interaction V = —It. First, 
the right panels of both figures reveal that the characteristic shapes of the boundary 
LDOS D{j = 1, e) in figures [6] and [7] are purely local features. The LDOS D{j > 2, e) on 
other sites look completely different, even on the closest sites. Therefore, it is unlikely 
that the boundary exponent (fT^ can be determine precisely in STS experiments which 
average over several A [TT]. Second, we note that the LDOS for j > 1 are quite similar 
for attractive and repulsive interactions besides the different band width ^ At*. The 
main reason for this likeness is that the spectra D{j, e) are clearly dominate by j peaks 
up to (at least) j = 10 for all interactions V. These oscillations are similar to those 
found in the LDOS for non-interacting spinless fermions (fT5|) . Finally, in the left panels 
of figures M and [9] we see that the amplitude of these oscillations decrease rapidly with 
increasing distance j from the boundary. Already for j ^ 50 = N/4 both LDOS do not 
appear to change on the scale of the broadening rj = 0.08t. This confirms that DDMRG 
calculations of the LDOS in the middle of a chain yield a (broadened) bulk density of 
states. 

The above discussion is based on our DDMRG data in chains with up to 200 sites 
and T] > O.OSt. As seen in the study of bulk and boundary density of states, this 
resolution is not good enough to analyze the asymptotic behaviour close to the Fermi 
energy. Moreover, the multiple oscillations observed in the LDOS in the crossover regime 
1 < j ^ imply that a power-law behaviour could only be observed in a significantly 
smaller energy range than for the bulk and boundary density of states. We estimate 
that carrying out such a study would require a computational effort which is one order 
of magnitude higher than for the present work (which amounts to 2 ■ lO'^ CPU hours). 
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Figure 9. LDOS D{j, e) calculated on a 200-site chain (77 = O.OSt) for V = -t. (a) 
From one edge to the middle of the chain: j = 1, 10, 20, . . . , 100 from top to bottom, 
(b) Close to the chain edge: j = 1, 2, . . . , 10 from top to bottom. 



The crossover from boundary to bulk density of states in Luttinger liquids has 
already been studied with field theory |2UI2S]. We have not performed any quantitative 
comparison between DDMRG spectra and the field-theoretical results or the fRG 
results |26] because of the limited resolution of our data. Nevertheless, there is a clear 
discrepancy between our results and one of the main predictions of bosonization. In an 
extension of the Tomonaga-Luttinger model with hard wall boundaries it has been shown 
that boundaries affect the LDOS D{j, e) only below a crossover energy [21]. Above this 
energy bulk behaviour is observed. This crossover energy is given by vh/{ja) = 2t*/j. 
However, in DDMRG spectra (see figures M and [9]) boundary effects are clearly seen as 
rapid oscillations over the full band width, even for j = 20. Thus our results indicate 
that there is no crossover energy above which only bulk behaviour is observed in the 
spinless fermion model. This is also demonstrated by the exact LDOS of noninteracting 
fermions (|T5|) which can be regarded as a Luttinger liquid with K = 1. A correct 
statement for the spinless fermion model is that boundary effects are only visible in the 
LDOS D{j,e) at a resolution < 4:t*/j. If the LDOS is examined at a resolution larger 
than this value (for instance, if it has been convolved with a Lorentzian or Gaussian 
distribution of width > 4t*/j), then boundary effects are not visible. Thus in the 
DDMRG spectra boundary effects disappear when the broadening (flOl) exceeds a value 

PS 4r/j. 

5. Conclusion 

We have investigated the local density of states (LDOS) -D(j, e) of the one- dimensional 
half-filled spinless fermion model in its Luttinger liquid phase. Using the DDMRG 
method we have determined the LDOS over its full energy range with a resolution of 
7] = O.OSt or better. In all cases most of the spectral weight is concentrated in the band 
|e| < 2t* of elementary excitations from the Bethe Ansatz solution. We have found that 
the bulk density of states can be identified with the LDOS computed in the middle 
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of a chain if we use a large enough broadening rj. It is dominated by sharp peaks at 
the band edges e = ±2t* for all interactions V. For repulsive interactions {V > 0) the 
boundary LDOS (on the first lattice site) consists in two symmetric joined bands which 
are reminiscent of Hubbard bands. For attractive interactions {V < 0) the boundary 
LDOS is dominated by a peak at the Fermi energy, which becomes increasingly sharp 
as V decreases and turns into a Dirac (5-peak for V — j- —2t. 

We have also examined the bulk and boundary density of states close to the Fermi 
energy with a resolution rj = O.Olt on lattices with 1600 sites. In the bulk density of 
states for V ^ and in the boundary LDOS for V > we can observe the power-law 
pseudogap D{j,e) ~ |e|" predicted by the Luttinger liquid theory over an energy range 
~ t* although we have to compensate explicitly for finite-size effects when the exponent 
a is large. In the boundary LDOS for < we see the power-law divergence predicted 
by the Luttinger liquid theory. Using a finite-size scaling analysis of DDMRG data on 
lattices with up to 3200 sites we obtain exponents a which agree quantitatively with 
the exact results given by the Luttinger liquid theory combined with the Bethe Ansatz 
solution. 

Finally, we have explored the crossover from the boundary LDOS D{j = l,e) to 
the bulk density of states D{j = N/2,e) as one moves away from the chain edge. The 
characteristic features of the boundary LDOS D{j = 1, e) in Luttinger liquids appear to 
be a purely local effect as all D{j > 2, e) look completely different. They are dominated 
by oscillations which are similar to those observed in the non-interacting chain. The 
crossover is smooth and boundary effects can be visible far away from the chain edge 
(j ^ 1) at all energies e but vanish beyond a point set by the DDMRG resolution, i.e. 
when 4t*/j < rj. 

The LDOS of the spinless fermion model cannot be used directly to interpret 
photoemission or STS experiments. For this purpose we have to consider electronic 
models such as the Hubbard model. The present work has been motivated in part by 
the failure to observe a power law in the density of states of the Hubbard model using 
DDMRG with an energy resolution rj = O.lt. For spinless fermion models the DDMRG 
method allows us to study larger system sizes and thus to reach a better resolution 
~ 1/A^ than for electronic models. Our work confirms that the power-law features of 
Luttinger liquid LDOS can be studied with DDMRG. However, the resolution required 
for the Hubbard model is certainly much smaller than O.lt. In the spinless fermion 
model we have to use a resolution of about lO'^t* for a ^ 1 and lO^^t* for a ^ | to 
observe a power law. In the Hubbard model a < 1/8 and the energy scale of the power- 
law behaviour is much lower than the bare band width 4t. This scale is exponentially 
small ~ exp{—t/U) for small U/t [18] and it is probably set by the effective exchange 
coupling J ~ t'^/U for large U/t. Thus a resolution even smaller than O.Olt is required 
which implies system sizes well in excess of 10'^ sites. 

The recent experimental observation of Luttinger liquid LDOS in atomic wires 
certainly calls for a renewed effort in the study of the LDOS in Hubbard-type models 
for interacting electrons in one dimension. Fortunately, larger exponents a are achieved 
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in extended Hubbard models with a non-local Coulomb interaction [271 EH], which 
are anyway more realistic than the original Hubbard model. We think that in these 
electronic models the LDOS and the power-law pseudogap could be investigated with a 
sufficient accuracy and a reasonable computational effort using the DDMRG method. 
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